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SUMMARY 

A computational methodology is described to probabilistically simulate the stress concentra- 
tion factors (SCF’s) in composite laminates. This new approach consists of coupling probabilistic 
composite mechanics with probabilistic finite element structural analysis. The composite mechanics 
is used to probabilistically describe all the uncertainties inherent in composite material properties, 
whereas the finite element is used to probabilistically describe the uncertainties associated with meth- 
ods to experimentally evaluate SCF’s, such as loads, geometry, and supports. The effectiveness of the 
methodology is demonstrated by using it to simulate the SCF’s in three different composite laminates. 
Simulated results match experimental data for probability density and for cumulative distribution func- 
tions. The sensitivity factors indicate that the SCF’s are influenced by local stiffness variables, by load 
eccentricities, and by initial stress fields. 


INTRODUCTION 

It is generally accepted that flawed structures (those with holes or cracks) fail because stress 
concentrations cause damage of such magnitude that (1) the structure cannot safely perform as 
designed and qualified or (2) catastrophic global fracture is imminent. According to the authors’ 
knowledge, this is true for structures made from traditional homogeneous materials as well as for fiber 
composites. The difference between fiber composites and traditional materials is that composites have 
multiple fracture modes that initiate local flaws whereas traditional materials have only a few. Stress 
concentrations in composites are influenced by uncertainties in many more factors than in metals 
because of the inherent anisotropic and layered structure of composites. Any predictive approach for 
simulating stress concentrations in fiber composites needs to formally quantify (1) all possible fracture 
modes, (2) the types of flaws they initiate, and (3) the uncertainties that influence the magnitude of the 
stress concentration. 

One of the ongoing research activities at the NASA Lewis Research Center is the development 
of a methodology for “probabilistic structural analysis” in general. A part of this methodology consists 
of step-by-step procedures to probabilistically quantify the uncertainties in composite behavior such as 
ply and laminate mechanical, thermal, and heat transfer properties (refs. 1 and 2) and composite 
structural response such as displacement, stress, and vibration frequency (ref. 3). The methodology has 
been embedded in the computer code IPACS (Integrated Probabilistic Assessment of Composite 
Structures, ref. 4) through the integration of two NASA Lewis in-house codes: PICAN (Probabilistic 
Integrated Composite ANalyzer, ref. 5) and NESSUS (Numerical Evaluation of Stochastic Structures 
under Stress, ref. 6). Another new application of that approach is to probabilistically simulate the range 



of scatter (uncertainties) in stress concentration factors (SCF’s) in composite laminates. This simulation 
provides a formal method to (1) assess alternative designs with respect to their reliability and/or (2) to 
quantify the probability that the magnitude of the SCF’s would be within an acceptable range. This 
report describes the fundamental aspects of this approach and illustrates its application by probabilis- 
tically simulating SCF’s in several composite laminates and discussing their significance. 


SYMBOLS 


e eccentricity 

M moisture 

P pressure 

T use temperature 

t time 

X i design parameter 

O stress 

O a applied stress 

Gq stress in composite 

G x axial stress 

shear stress 

a y transverse stress 


FUNDAMENTAL CONCEPTS 
Computational Simulation 

A brief description of the fundamental concepts for the computational methodology is included 
for completeness. Uncertainties (scatter) in composite laminates and composite structures result from 
their inherent nature and from their fabrication process (fig. 1). The SCF’s in composite laminates are 
probabilistically assessed by using the IP ACS code (ref. 4). This code integrates the probabilistic 
composite mechanics in PICAN (ref. 5) with the probabilistic finite element structural analysis in 
NESSUS (ref. 6). A schematic of this integrated computer code is shown in figure 2. 

Probabilistic composite mechanics is used to formally describe all the uncertainties associated 
with the composite — from micromechanics to laminate. The uncertainties are incorporated in the 
primitive variables that describe the laminate in the composite mechanics equations. These include 
constituent material properties, fabrication process variables, ply orientation angles, and ply thickness 
(refs. 1 and 2). These uncertainties are accounted for in PICAN (fig. 2). Typical values for boron, 
glass, graphite, and epoxy materials are listed in table I. 

Probabilistic structural analysis is used to formally describe the uncertainties associated with 
the structure. The uncertainties are incorporated in the primitive variables in the structural mechanics 
equations that define the structure. They include special geometry, boundary conditions, and load 
conditions (ref. 3). These uncertainties are accounted for by NESSUS (fig. 2). 

The simulation process in IPACS is as follows: 

(1) The scatter in all the primitive variables that describe the composite can be represented by 
well-known probabilistic distributions. 

(2) The values for the primitive variables can be randomly selected from the known distribu- 
tions for a specific composite. 
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(3) These values can be used in composite mechanics to predict composite behavior. 

(4) The whole process can be repeated many times to obtain sufficient information to develop 
the distribution of the ply properties, composite properties, or structural responses. 

These concepts are analogous to making and testing composites. The probabilistic distributions 
represent uncertainties of available materials that the composite can be made from. The composite 
mechanics represents the physical experiment, and the process repetition represents several experi- 
ments. Subsequent statistical analysis of the data is the same for both approaches. The primitive varia- 
bles that describe the composite are identified by examining the fabrication process, the structural 
assembly, the boundary conditions, and the loading conditions. During the simulation process, sensitiv- 
ity factors are also evaluated. The sensitivity factors rank all the participating random variables 
according to the probable significance of a specific structural response. The commonly used sensitivity 
in a deterministic analysis is the performance sensitivity dZ/dX^ which measures the change in the 
performance Z due to the change in a design parameter X-. This concept is extended to the probabilis- 
tic simulation to define the probabilistic sensitivity, which measures the change in the probability or 
reliability relative to the change in each random variable, including both physical significance and 
scatter of that random (primitive) variable. 


Closed-Form Equation 

Another fun dam ental concept includes the closed-form equations (CFE s) for concentration 
factors (ref. 7). The equations are based on circular holes in infinite laminates. The equations for the 
SCF’s are included in PICAN, which can be used to evaluate them probabilistically. Then, the results 
can be compared with those obtained from IP ACS. 

RESULTS AND DISCUSSION— LAMINATE STRESS CONCENTRATION FACTORS 

SCF’s for three different laminate composite systems (boron/epoxy. E-glass/epoxy, and 
graphite/epoxy, respectively) were evaluated. The finite element model used, the various uncertainties 
considered (including combined stresses, results obtained, and their significance to laminate integrity), 
and comparisons with the CFE and with available data are discussed in this section. 


Laminate Finite Element Model 

The laminates were 20 by 10 in. with a center hole of 1 in. diameter and a thickness of 
10 plies, or about 0.05 in. They were loaded (1) in tension parallel to the 20-in. dimension and 
(2) under combined stresses to evaluate these effects on the SCF. Figure 3 shows the finite element 
model of the laminat es. The laminate configuration, [0/±45/0/] s , was the same for all three composite 
systems. 


Boron Epoxy Laminate 

Figure 4 shows the scatter in the SCF in terms of the probability density function (PDF), or 
frequency of occurrence. The range of uncertainties included in the IPACS simulation are summa- 
rized in the figure for the reader’s convenience. Three different distributions are shown in this figure: 
(1) simulated by IPACS, (2) independent source (the Monte Carlo simulation), and (3) experimental 
data The independent source distribution and the experimental data are from reference 8. The PDF 
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and the cumulative distribution function (CDF) for the experimental data were generated by the 
authors from the range and mean value of the experimental data given in reference 8. IPACS simulates 
the experimental data almost exactly, whereas the independent source overpredicts the mean and the 
range. Figure 4 also shows the corresponding CDF (probability of occurrence). Again, the IPACS and 
experimental data are identical for all practical purposes. 

In addition, figure 4 shows the sensitivity factors for two levels of probability. It is interesting 
to note that of about 40 primitive variables used, only 4 have a significant effect on the SCF. All four 
of these primitive variables define local stiffness, demonstrating that SCF’s are predominantly stiffness 
controlled. Comparisons of the IPACS’ simulated range of uncertainty in the SCF with that predicted 
by using the CFE (ref. 7) are shown in figure 5 for two ranges of uncertainties in the primitive varia- 
bles. Note that the experimental da t a range is shown by the cross-hatched bar. The following observa- 
tions are worth noting: 

(1) The CFE overpredicts the mean of the stress concentration by about 20 percent. 

(2) The CFE predicts a smaller range for both the 2 and 5 percent uncertainties than the 
IPACS simulations do. 

(3) The CFE does not agree with the experimental data. 

(4) A 5-percent uncertainty in the primitive variables increases the uncertainty range by almost 
4 times in the IPACS simulation but only by 2 times in the CFE. 

An important conclusion is that the whole laminate must be modeled in order to adequately 
represent the uncertainty range in the SCF. It is very important to keep in mind that the uncertainty 
ranges simulated by the SCF’s must be compared with the laminate’s respective local strength in order 
to assess this laminate’s resistance to fracture initiation for a specific applied stress magnitude. 


Glass Epoxy Laminate 

Figure 6 shows the uncertainty range of the SCF in glass/epoxy laminates. Note that a 
3-percent uncertainty of the laminate thickness out-of-plane eccentricity was included in this laminate. 
The reason for this eccentricity is that the IPACS simulation predicted the shape of the PDF but not 
the mean of the experimental data (ref. 8). The mean without the eccentricity was about 3.2 compared 
with a mean of about 3.5 for the experimental data. Other uncertainties such as in-plane eccentricities 
and greater range in ply misorientation had negligible influence on the mean. Hindsight shows that 
probable eccentricity should have been included in the initial simulation. An important observation is 
that IPACS captures the uncertainty range of the SCF provided that all probable uncertainties are 
included. As a side note, the same amount of eccentricity had a negligible effect on the SCF for the 
boron/epoxy laminate. 

Figure 6 also shows the corresponding CDF for the glass/epoxy laminate. For all practical 
purposes, the IPACS simulation reproduces the experimental data. 

The sensitivity factors are summarized in the bar chart in figure 6. Note that the sensitivity 
factor for the out-of-plane eccentricity is of about the same order as that for the ply thickness. The 
eccentricity is fourth in significance after the longitudinal fiber modulus, the fiber volume ratio, and 
the ply misorientation range. Uncertainties in the remaining primitive variables for the glass, fiber, 
matrix and the fabrication variables (table I) have a negligible effect. 

An important conclusion from the collective results shown in figure 6 is that the probable 
uncertainties in all participating primitive variables need to be included in the probabilistic simulation 
of a specific laminate property or response. It is generally prudent to do so because the sensitivity 
factors will rank the significance of each primitive variable and the dominant ones will be identified. 
These dominant variables can then be used in subsequent detailed analysis as required. 
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Graphite Epoxy Laminate 


Figure 7 shows the uncertainty range in the SCF of graphite/epoxy laminates in comparison 
to experimental data and predictions from an independent source (ref. 8). As for the glass-epoxy 
laminates, a 7-percent out-of-plane eccentricity of the laminate thickness was used to shift the IPACS- 
simulated PDF to the right to match the experimental data. Otherwise, the shape of the PDF was 
identical. 

Figure 7 also shows the corresponding CDF’s. They are coincident as would be expected from 
the very good agreement of the PDF’s. The respective sensitivity factors are shown in the bar chart in 
this figure. The three dominant primitive variables are those that control local stiffness — fiber modulus, 
fiber volume ratio, and ply misorientation. The eccentricity is the fourth dominant primitive variable 
with a magnitude of about one-half of the fiber modulus. Two obvious questions surface at this point 
for both the glass/epoxy and the graphite/epoxy laminates: 

(1) If the eccentricity is not dominating (does not have a high sensitivity factor), why should it 

be even included? 

(2) Why not increase the range in the dominating variables? 

The answers follow: 

(1) The sensitivity factors of the dominant primitive variables include both physical and 
probabilistic effects. 

(2) Increasing the range of the dominating primitive variables spreads the uncertainty range of 
the SCF but does not shift the mean. 

These answers are important and lead to the following general observation: When the shapes 
of the IPACS-simulated probabilistic distributions are identical, or nearly so, with those of the mea- 
sured data but are separated by some amount, then either (1) the in situ mean of one of the primitive 
variables is off (by an amount equal to the difference between the two means) or (2) a primitive varia- 
ble is missing in the probabilistic simulation. This can be generally resolved by examining the specific 
case and may require additional simulations (such as a finer mesh near the hold if the probabilistically 
simulated mean is less than that of the data). The feedback from the results will indicate which is 
correct. 


Uncertainties of Initial Thermal Stress Fields 

It is well known that because of their fabrication process composite laminates have initial 
thermal stress fields that require cooling down from cure temperatures to room temperature during 
curing. It is reasonable to assume that this temperature will not be uniform and will cause nonuniform- 
ities in the residual stresses. These nonuniformities, in turn, may cause eccentricities and affect the 
laminate SCF. Although these stress nonuniformities are not known, they can be represented as 
uncertainties and quantified by using IP ACS to simulate their effects on the SCF. 

The effects on initial stresses due to a temperature through-the-thickness gradient of 80 °F with 
a ±5 °F scatter on the SCF in glass/epoxy laminates are shown in figure 8. All the other uncertainties 
were identical except the uncertainty for the out-of-plane eccentricity. As can be seen, IPACS 
simulates the uncertainty range in the SCF without requiring out-of-plane eccentricity. This is 
significant indeed because the uncertainties due to initial thermal stresses are more likely to be unsus- 
pected than those caused by load eccentricities. Comparing the sensitivity factors to those in figure 6 
shows that the dominance of the fiber modulus increased from 0.6 to 0.7, whereas the dominance of 
the fiber volume ratio decreased from 0.6 to 0.5. The temperature replaced the dominance of the 
eccentricity with approximately the same magnitude. This observation implies that the thermal stress 
effect causes la minafp out-of-plane displacements comparable to those induced by eccentricities. 
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Corresponding initial stress uncertainties in the SCF of graphite epoxy laminates are shown in 
figure 9. An important observation is that the initial stress uncertainties are not sufficient to eliminate 
the out-of-plane eccentricity for this laminate. However, these uncertainties reduced the eccentricity’s 
significance from 7 percent to 3.5 percent, or by one-half. The sensitivity factors show that the temper- 
ature uncertainty is of about the same magnitude as the out-of-plane eccentricity uncertainty. As for 
the glass/epoxy laminates, the modulus sensitivity increased in significance, whereas the fiber volume 
ratio sensitivity decreased. Another important conclusion is that temperature is another primitive varia- 
ble which introduces uncertainties that must be included in the IPACS simulation. 


Uncertainties of Combined Stresses 

The combined stresses that may result from the test fixture are another source of probable 
uncertainties. These effects are shown in figure 10 for the graphite/epoxy laminate. Some important 
observations are that the uncertainties of o y and stresses have no influence on the range of uncer- 
tainties in the SCF; however, the uncertainty of the a x stress changes the SCF’s shape. 

The corresponding sensitivity factors are also shown in figure 10. The sensitivity factor for the 
o x stress dominates and even exceeds that for the fiber modulus, which was previously the dominant 
primitive variable for the SCF. Also, the sensitivity factor for the o y stress is greater than those for 
the ply thickness and eccentricity. However, the residual stress due to cure temperature did not even 
appear in the ranking. One conclusion is that the residual stress affects only ply stresses that are inter- 
nal to the laminate and that cannot be evaluated from laminate tests. Another observation is that the 
applied stress field effect, depending on its magnitude, may override the effect of any eccentricities 
from residual stresses. Back-to-back strain gauges will distinguish this effect in the experimental eval- 
uations especially under low levels of load. 


GENERAL COMMENTS 

From the previous observations and discussion, the following general comments are worthy of 
note for the probabilistic assessment of SCF. 

(1) Probabilistic evaluation provides a priori information for preparing a test plan and esti- 
mating the range of measured values of the SCF if the laminate was (a) fabricated as designed and 
(b) tested according to the test plan. 

(2) Sensitivity factors offer a direct means to identify which primitive variables dominate, 
which variables are insignificant, and which may have been missed in the simulation. They also 
provide guidelines to improve designs and may indicate what went wrong during testing. 

(3) Probabilistic evaluation provides guidelines for (a) the minimum number of experiments 
and (b) quality criteria for material acceptance and for fabrication tolerances to assure that specific 
designs will meet reliability and safety requirements. 

(4) Probabilistic evaluations combined with probabilistic strengths will provide quantifiable 
information to assess risk for a specified reliability. 

Collectively, these points provide designers and manufacturers with metrics to make judicious 
decisions that result in reliable products with specified life-cycle costs. 


SUMMARY OF RESULTS 

The stress concentration factor (SCF) in composite laminates was probabilistically evaluated 
for three different composite systems (boron, glass, and graphite fiber (epoxy). The SCF’s were 


6 



simulated computionaUy with the computer code IPACS (Integrated Probabilistic Assessment of 
Composite Structures), which consists of probabilistic composite mechanics and probabilistic finite 
element structural analysis. The following results were obtained: 

1. The SCF means and scatter ranges predicted by IPACS are in very good agreement with 
available data for the three laminates. 

2. The sensitivity factors indicate that the primitive variables for local stiffness (fiber, modulus, 
fiber volume ratio, and ply misorientation) have the greatest influence on the probable magnitude of 
SCFs. 

3 a 7 -percent out-of-plane eccentricity in the load application increased the mean by about 
22 percent of the SCF in graphite/epoxy laminates. Back-to-back strain gages must be used to measure 
SCF’s in composite laminates. 

4. Nonuniformities in initial stress fields (for example, unsymmetric residual stress induced by 
curing) cause laminate bending. The effects of these nonuniformities on SCF’s are similar to those of 
out-of-plane load eccentricities. 
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TABLE L — PRIMITIVE VARIABLES AND ASSOCIATED UNCERTAINTIES 

(a) Fibers. 


Primitive variable 

Units 

Distribution 

type 

Mean value 

Fiber 

Boron* 

Glass b 

Graphite b 

Modulus 








In longitudinal direction 

mpsi 

Normal 

58.0 

12.4 

38.0 

In transverse direction 





58.0 

12.4 

1.1 

Shear modulus 








In plane 





24.2 

5.17 

1.5 

Out of plane 

y 




24.2 

5.17 

.8 

Poisson’s ratio 








In plane 

inVin. 



0.2 

.2 

.2 

Out of plane 

inVin. 



0.2 

.2 

.25 

Thermal expansion coefficient 








In longitudinal direction 

ppm/°F 



2.8 

2.8 

-.55 

In transverse direction 

ppm/°F 



2.8 

2.8 

5.6 

Density 

lb/in. 3 



.095 

.09 

.063 

Number of fibers per end 

- 


Fixed 

1 

204 

10 000 

Diameter 

in. 

Normal 

.004 

.00036 

.0003 

Heat capacity 

BTU/lb/°F 



.31 

.17 

.17 

Heat conductivity 








In longitudinal direction 

BTU-in./hr/in//°F 



1.55 

.625 

4.03 

In transverse direction 

BTUin./hr/in. 2 /°F 



1.55 

.625 

.403 

In out-of-plane direction 

BTU -inThi/in. 2 /°F 



1.55 

.625 

.403 

Strength 








Tensile 

ksi 

Weibull 

600 

360 

350 

Compressive 

ksi 

Weibull 

700 

60 

250 


(b) Matrix and fabrication variables. 


Primitive variable 


Modulus 
Shear modulus 
Poisson’s ratio 

Thermal expansion coefficient 

Density 

Heat capacity 

Heat conductivity 

Tensile strength 

Compressive strength 

Shear strength 

Moisture coefficient 

Diffusivity 


Fiber volume ratio 
Void volume ratio 
Ply thickness 
Ply misalignment 


Units 

Distribution 


type 


Mean value 


Matrix* 


mpsi 

Normal 

0.5 

mpsi 



.185 

inVin. 



.35 

ppm/°F 



42.8 

lb/in. 3 



.0443 

BTU/lb/°F 



.25 

BTUin7hriin. 2 /°F 



.104 

ksi 

Weibull 

15.0 

ksi 

Weibull 

35.0 

ksi 

Weibull 

13.0 

in7in71 -percent 

Normal 

.004 

moisture 

Normal 

.0002 


Fabrication variables 


percent 

Normal 

*50 

percent 



*0.0 

in. 



*.0055 

in. 



c 0 
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Fabrication schematic 

• Constituents 

(Fiber and matrix properties) 

• Fiber misalignment 

• Fiber volume ratio 

• Void volume ratio 

• Ply misalignment 

• Ply thickness 

Figure 1 . — Random variables in the typical fabrication 
process. 
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Figure 3. — Finite element modeling. 
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IPACS 

Monte Carlo simulation (ref. 8) 

Experimental data (ref. 8) 



VV= 10 in. 


SCF, oa / 



■■ Fiber modulus 
rm Fiber volume ratio 
ES Ply misalignment 
F3S3 Ply thickness 




At 0.001 probability At 0.999 probability 


SCF, cr^/ac 


Figure 4.— Probabilistic stress concentration factor (SCF) for a (0/±45/0/90)s boron/epoxy laminate plate. Fiber and matrix properties with 
2.5-percent scatter; thickness with 2.5-percent scatter; and ply alignment with 2.25° scatter. 




SCF 

Figure 5. — Probabilistic stress concentration factor (SCF) for a 
(0/±45/0/90)s boron/epoxy laminate plate. 
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Cumulative distribution function, CDF 


IPACS 

Monte Carlo simulation (ref. 8) 

Experimental data (ref. 8) 



SCF, cr^/crc 

Figure 6.— Probabilistic stress concentration factor (SCF) for a (0/±45/0/90)s glass/epoxy laminate plate. Fiber and matrix properties with 
1 -percent scatter; thickness with 1 -percent scatter; ply alignment with 0.9° scatter; mean of eccentricity, e, 3 percent of Th, and 
coefficient of variation of e, 2 percent of mean. 


Cumulative distribution function, CDF 


C 



3.5 4.0 4.5 At 0.001 probability At 0.999 probability 

SCF, c 

Figure 7.— Probabilistic stress concentration factor (SCF) for a (0/±45/0/90)s graphite/epoxy laminate plate. Fiber and matrix properties 
with 1 -percent scatter; thickness with 1 -percent scatter; ply alignment with 0.9° scatter; mean of e, 7 percent of Th; and coefficient of 
variation of e, 2 percent of mean. 
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Cumulative distribution function, CDF 



SCF, OfljGQ 


igure 8.— Probabilistic stress concentration factor (SCF) for a (0/±45/0/90)s glass/epoxy laminate plate (initial stress effect). Fiber and 
matrix properties with 1 -percent scatter; thickness with 1 -percent scatter; ply alignment with 0.9° scatter; and reference temperature 
70 °F, and cure temperature, 150 °F, with 5 °F scatter. 



Cumulative distribution function, CDF 



3.5 


4.0 

SCF, c JA^C 


4.5 


At 0.001 probability 


At 0.999 probability 


Figure 9. — Probabilistic stress concentration factor (SCF) for a (0/±45/0/90)s graphite/epoxy laminate plate (initial stress effect). Fiber and 
matrix properties with 1 -percent scatter; thickness with t -percent scatter, ply alignment with 0.9° scatter; reference temperature, 70 °F, 
and cure temperature, 1 50 °F, with 5 °F scatter; mean of e, 3.5 percent of Th; and coefficient of variation of e, 2 percent of mean. 
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Cumulative distribution function, CDF 





SCF, o'^/ctc 

Fiqure 10.— Probabilistic stress concentration factor (SCF) for a (0/±45/0/90)s graphite/epoxy laminate plate (combined stresses effect). 
Fiber and matrix properties with 1 -percent scatter; thickness with 1 -percent scatter; ply alignment with 0.9° scatter; reference 
temperature, 70 °F, and cure temperature, 350 °F, with 10 °F scatter e, 3.5 percent of 7 h, with 2-percent scatter and <r x , <x y , and 

with 1 -percent scatter. 


17 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden lor this collection erf information is estimated to average 1 hour per response, including the time tor reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway. Suite 1204, Arlington. VA 22202-4302. and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 


4. TITLE AND SUBTTTLE 


2. REPORT DATE 


December 1994 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Probabilistic Simulation of Stress Concentration in Composite Laminates 


6. AUTHOR(S) 


C.C Chamis, P.L.N. Murthy, and D.G. Liaw 


WU-5 10-02-12 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-9077 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 1067 12 


11. SUPPLEMENTARY NOTES 


Prepared for the 34th Structures, Structural Dynamics, and Materials (SDM) Conference cosponsored by AIAA, ASME, 
ASCE, AHS, and ASC, La Jolla, California, April 19-22, 1993. C.C. Chamis and P.L.N. Murthy, NASA Lewis Research 
Center; D.G. Liaw, Sverdrup Technology, Inc., Lewis Research Center Group, Brook Park, Ohio 44142 (work funded by 
NASA Contract NAS3-25266). Responsible person, C.C. Chamis, organization code 5200, (216) 433-3252. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 39 


13. ABSTRACT (Maximum 200 words) 

A computational methodology is described to probabilistically simulate the stress concentration factors (SCFs) 
in composite laminates. This new approach consists of coupling probabilistic composite mechanics with probabilistic 
finite element structural analysis. The composite mechanics is used to probabilistically describe all the uncertainties 
inherent in composite material properties, whereas the finite element is used to probabilistically describe the uncertainties 
associated with methods to experimentally evaluate SCFs, such as loads, geometry, and supports. The effectiveness of 
the methodology is demonstrated by using it to simulate the SCFs in three different composite laminates. Simulated 
results match experimental data for probability density and for cumulative distribution functions. The sensitivity factors 
indicate that the SCFs are influenced by local stiffness variables, by load eccentricities, and by initial stress fields. 



14. SUBJECT TERMS 

Probabilistic mechanics; Probabilistic finite-element; Constituent materials; Fibers; 
Matrices; Probability distributions; Eccentricities; Thermal stress; Combined stresses; 
Sensitivities; Boron; Glass; Graphite; Epoxy 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


15. NUMBER OF PAGES 

19 


16. PRICE CODE 


20. LIMITATION OF ABSTRACT 


NSN 7540*01 -280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 




















Jg 
is 
-? e - 
- 03 

o C 
^ at 
•V 5 ' 
9- (D 
< at 
a> at 


8 

8 


o 

33 


© 

c 

3 


National Aeronautics and 
Space Administration 

Lewis Research Center 

21000 Brookpark Rd. 
Cleveland, OH 44135-3191 



